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Abstract 

Renormalization group and Coulomb gas mappings are used to derive theoretical 
predictions for the corrections to the exactly known asymptotic fractal masses of 
the hull, external perimeter, singly connected bonds and total mass of the Fortuin- 
Kasteleyn clusters for two-dimensional (/-state Potts models at criticality. For q = 4 
these include exact logarithmic (as well as log log) corrections. 



1. INTRODUCTION 

(/-state Potts models, with interaction —J5 (7i<Tj (cjjj = 1,2, ...,q) for the nearest neighbor 
(nn) sites i,j, have played an important role in condensed matter physics^. Here we study 
geometrical aspects of the critical Potts clusters, in two dimensions. For an arbitrary con- 
figuration of Potts states, one creates bonds between neighboring sites which have the same 



l 
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state, Oi = <jj, with a probability p = 1 — exp (—J/kT). No bonds are created between sites 
with <Tj ^ <7j. Here we study the fractal geometry at T c of the clusters, made of sites con- 
nected by bonds!. Specifically, we measure the fractal dimensions Dm, Djj, Drp, and D$c 
describing the scaling of the cluster's mass, hull, external accessible perimeterB and singly 
connected bonds!, respectively, with its radius of gyration R. As emphasized by Coniglicd, 
many of these fractal dimensions have been derived analytically!. Some others have been 
found more recently!'!. These theoretical values are summarized in Table 1 , in terms of the 
Coulomb gas coupling constant g = ^ arccos(— ^). 



Table 1: Exact theoretical predictions. 
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We are currently studying the geometry of such Potts critical clusters, using numerical 
Monte Carlo simulationaHI. These simulations show that the asymptotic power law depen- 
dence of the various masses on R is approached relatively slowly, and therefore the analysis 
of the data must include correction terms, particularly as q approaches 4. The present paper 
contains a brief summary of analytic results for several of these corrections. Our predictions 
for the correction exponents (apart from the analytic ones) are also listed in Table 1. 



2. RENORMALIZATION GROUP 

The first correction relates to the dilution field ip, which is generated under renormalization 
even when one starts with the non-dilute casdM For our non-dilute case, one expects ip(£) to 
increase under the renormalization group recursion relations (RGRR's) from — oo towards 
its critical value ip* (< for q < 4), while the cluster linear size, like all other lengths, 
rescales as R(£) = Re~^. Following Cardy et a/.0, we assume that t§ prefacing iterations 
bring tp from — oo up to VK^o) = "00 > with |-0o| *C 1. We then expand the RGRR for ip in 
powers of ifi and e = q — 4, 

= a(e + ip 2 + bip 3 + reip + ...). (1) 

At q = 4, this yields 

a{i - 4) « 1 " + blog b + m S l 1 ■ (2) 
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Iterating up to R(£ x ) = 1, we can express ip(£ x ) in terms of logi? = £ x . For large R, 
V'o/W*) ~ A(logR+Blog(logR) + E) + 0(loglogR/logR), where A = -atp , B = -b/a 
and E also depends on a, 6, £q and Tpo. 

For q < 4, an expansion to second order in e' = \/— e yields ip(£ x ) « ■0* + BR~ e , with 
0* = - e '(l - ( r - 6)e'/2) + . . ., 9 = 2ae'(l - be') + . . . and B oc (0o - V*)- To order-e', one 
obtains a full solution, 

1 + B e -oi 

where B = (ip + e')/(0o ~~ O- Indeed, ip approaches ip* for large £. 

To obtain the scaling of Mg(R), we write the RGRR for the field hs conjugate to the 
density ps = Ms/R d as 

^f = (ys + c 5^(l + e s V + /s^ 2 + • • -))h s , (4) 

where the coefficients may depend on e. ps is then found as a derivative of the free energy 
with respect to hs- For q = 4, its singular part becomes 

Ps(^) oc e- de h s (£)/h (5) 

= exp[(y 5 -d)£+ / (c s ^(l + e s ^ + fs^ 2 + • • -))d£] 

oc e ( ys -^[^(£)/^o] Cs/a (l + 0(ip(l))). 

For large log R = £ x , this becomes 

M 5 oc fl Ds (log i? + Slog (log i?) + £)" C5/a (l + 0(loglo gj R/logi?)), (6) 

with Dg = ys(q = 4), and cs/a is to be taken from Table 1 (see below) . Note that 
B = —b/a is universal (i. e. independent of ipo), and the non-universal constant E is the 
same for all S. Equation (@) generalizes the logarithmic corrections of Cardy et aZ.0. 

In practice, the numerical results are always analyzed by looking at the local logarithmic 
slope, 

Df = dlogMs/dlogR = dlogh s /d£\ e=e x 

= ys + cs^(n{l + es^{e) + fs^) 2 + •••)• (7) 

In some cases, this expression (in which 0(£ x ) is related to logi? = £ x via Eq. (||)) gave a 
better fit than the derivative of the approximate expression in Eq. @. 
For q < 4, to leading order in e' , the same procedure turns Eq. (|3|) into 

M s oc R Ds (l - BR- d y cs / a R Ds (l + f s R~ 9 ), (8) 

where Ds ~ — c^e' and # ~ 2ae'. The RHS of this equation remains correct also for 
higher orders in e' . Note that to the lowest order in e', the ratios fs/fs' are universal, being 
equal to cs/cs>- This is similar to analogous ratios for thermodynamic properties in the 
usual e-expansion0. Expanding the exact Ds (Table 1) in e' yields c$. Using also a = 1/tt 
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(see below) yields our predictions for cs/a (given in Table 1), to be used in fitting Eq. (||). 
The form on the RHS of Eq. (|8|) is already implied by den NijH, who found that the 
pair correlation functions Gf/(r) can be expanded as a sum over r~ 2Xn , implying a leading 
correction exponent = 2(x n+ i — x ri ) = 4(4 — g)/g. Expanding this expression in powers of 
e' yields the coefficients a = 1/tt and b = — 1/2tt, which we use in our fits to Eq. (|2|). The 
value a = l/ir also agrees with Cardy et a/.0. This expression for 8 also reproduces known 
results for q = 2,3, as listed in Table 1. 



3. COULOMB GAS 

The second source of corrections involves new contributions to the relevant pair correlation 
functions in the Coulomb gas representations!^!. In some of the exact derivations, the (/-state 
Potts model renormalizes onto the vacuum phase of the Coulomb gas, involving 'particles' 
with electric and magnetic 'charges' (e, m). At criticality, the corresponding Coulomb gas 
has a basic 'charge' eft = \2—g/2\ mod 4. Various Potts model two-point correlation functions 
Gg(r) are then mapped onto Coulomb gas analogs, which give the probability of finding 
two charged particles at a distance r apart. Asymptotically, these are given by 



where 



Gu G w Nifr) oc r 2x l(n^M^m 2 )] (g) 
[(ei,mi),(e2,m2)J v ' ' V / 

CG _ e i e 2 gmim 2 

X [(ei,mi),(e 2 ,m 2 )] ~ ^ ^ ■ U U J 

Hence one identifies Dq = d — x9, G n , with d = 2. The results in Table 1 for 

° Lte ,mi),(e2,m2)J ' 

S = M were obtained by den Nijal3, who noted that the spin-spin correlation function of 
the Potts model maps onto a Coulomb gas total electric charge Q = —2(f), which splits into 
the two charges e\ t 2 = ±1 — 4> (and mi t 2 = 0). Continuing along similar routes, Saleur 
and Duplantieri used a mapping onto the body-centered solid-on-solid model, requiring a 
vort ex- ant i vortex pair with = —<j> and m^2 = ±1/2 or ±1 for the fractal dimensions of 
S = H or S = SC. The Table also contains Duplantier's recent result!! for Dep = 2 — x^p, 
which has not been expressed in terms of Coulomb charges. The results for x§ and x^(j are 
special cases of the expression xe = g£ 2 /32 — (4 — g) 2 /(2g), with 1 = 2 and 4 respectively^. 
For percolation (q = 1 and g = 8/3), this expression also yields x^p = ^3 = 2/3 for the 
external perimeter and Xq = xq = 35/12 for the gates to fjords.§ 

We now turn to corrections to the leading behavior, den Nij Jll derived such corrections for 
the order parameter correlation function. In that case, he noted that the charge Q = —2<p 
could also split into the pair ei^ = ±3 — <fi, yielding a contribution to Gfj of the form 

r 2x m,2^ w ith Xfr I2 = x f(3—^>o) (-3_0o)l = X M + V#- Since D = d — x usually represents a 
fractal dimension, we relate each of these correction terms to some subset of the cluster, 
with dimension Dm,2 = 2 — x^, I2 = Dm — 4/g. Writing Mm as a sum of powers i^^Jll we 
have M M oc R d m{1 + f'R- 6 '), with 6' = A/g. 

As far as we know, there has been no discussion of the analogous corrections to the other 
subsets discussed here. In the spirit of den NijsS, we note that the correlation function 
for both H and SC could also result from electrical charges e\ % = ±2 — <f>, instead of 
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—(f). For both of these cases this would give x' = x + 2/g, hence a correction exponent 
9" = 2/g. At the moment, there exists no theory for corrections to Mep- However, in 
the spirit of the renormalization group it is also reasonable to interpret 9' and 9" as the 
scaling exponent of some irrelevant perturbation (yet to be identified). If that were true 
then we might expect the same perturbation also to affect other quantities, like Mep- This 
conjecture is supported by the 'super universal' relation, (Dpr — 1)(Dep — 1) = 1/4, found 
by Duplantieri. If this relation also holds for the effective dimensions (as happens e. g. in 
the e-expansio n0), then H and EP should have the same correction exponents. 



The last source of corrections involves 'analytic' terms, coming e. g. from linear cuts with 
dimensions (Dg — 1), El or from replacing R by (R + a), since there are many possible 
candidates for the correct linear measure of the cluster. These would imply corrections of 
relative size 1/R. 

Combining all of these sources, we end up with the prediction (for q < 4) 



with 9{ = 9, 9' (or 9") and 1. Indeed, our numerical simulations^] basically confirm these 
expressions. 

In summary, we have presented several general expressions for the (/-dependent corrections 
to the asymptotic -R-dependence of the mass, hull, external perimeter and singly connected 
bonds. Such corrections are crucial for fitting numerical data. It would be nice to have a 
unifying theory, which would confirm these expressions in a rigorous way. It would also be 
nice to obtain similar corrections for other geometrical and physical quantities, e. g. the 
number of gates to fjordsi. 
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4. ANALYTIC CORRECTIONS; SUMMARY 



Df = D s + Y,fiR- 



(11) 
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